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ABSTRACT 

The properties of the short, energetic bursts recently observed from the 7-ray binary LS I +61°303, 
are typical of those showed by high magnetic field neutron stars, and thus provide a strong indication 
in favor of a neutron star being the compact object in the system. Here, we discuss the transitions 
among the states accessible to a neutron star in a system like LS I +61°303, such as the ejector, 
propeller and accretor phases, depending on the NS spin period, magnetic field and rate of mass 
captured. We show how the observed bolometric luminosity ( £ few x 10 35 erg s _1 ), and its broad- 
band spectral distribution, indicate that the compact object is most probably close to the transition 
between working as an ejector all along its orbit, and being powered by the propeller effect when it 
is close to the orbit periastron, in a so-called flip-flop state. By assessing the torques acting onto the 
compact object in the various states, we follow the spin evolution of the system, evaluating the time 
spent by the system in each of them. Even taking into account the constraint set by the observed 7-ray 
luminosity, wc found that the total age of the system is compatible with being w 5-10 kyr, comparable 
to the typical spin-down ages of high-field neutron stars. The results obtained are discussed in the 
context of the various evolutionary stages expected for a neutron star with a high mass companion. 
Subject headings: stars: magnetars - X-rays: binaries - X-rays individual (LS I +61°303) 



1. INTRODUCTION 

LS I +61°303 is one the few high-mass X-ray bi- 
naries (HMXB) discovered so far to e mit the largest 
part o f their luminosity at high energies (iHermsen et all 
fl977t IGregorv fc fa^iol \Wm lAlbert et a.1.1 I200(tD . be- 
ing therefore a member of the class of 7-ray bina- 
ries. Variability of its emission, at the timescale set by 
the ~26.5 d orbital p eriod, has been fou n d at almost 
all wa v elengths, e.g . , lAlbert et all (12008T) : lAbdo et al.l 
(pOOl : iTorres et alJ (|2010l ): iZhang et all ((20101) The 
companion star is a massive BOVe star, with a mass 
between 10 and 15 M Q , in an eccentric 26.5 d orbit 
(|Casares et al.l 120051) . For the nature of the compact 
object in 7-ray binaries, models involving an accret- 
ing black hole launching a relativistic jet (the micro- 
quasa r scenario; see, e.g. iBosch-Ramon fe Khangulvanl 
120091 and references therein) and a rotation-powered 
neutron star (NS in the follo wing) emitting a relativis- 
tic wind of particles (see, e.g. iMaraschi fe Treves! H98U: 
lDubusl[2006T iSierpowska-Bartosik fe Torresl 120081) . have 
been proposed. 

The presence of a NS in LS I +61°303 would be defi- 
nitely proven by the o bservation of pulsation s, but deep 
search es in the radio (McSwain et al.l l201lb and X-ray 
band (|Rea et al.ll2010b| ) were not successful, so far. This 
is not surprising, since free-free absorption easily washes 
out the pulses in the radio band, while the upper limit 
of « 10% (3 a confidence level) on the pulsed fraction 
in X-rays could well be larger than the actual pulsed 
fraction of the source. However, in the past few years, 
a couple of energetic (« 10 37 erg s _1 ), short ( £ 0.3 s) 
bursts were detected by the Swiff-Burst Alert Telescope 
from a region of a few arc minutes of radius, compatible 
with the position of LS I +61°303 (IDe Pasauale et al.1 
l2lMlBartbelmv et al.l[200llBurrows et al.ll2012L see Ta- 



blc[T]for their observed properties). The properties of the 
two bursts are typical of those observed in magnetars, 
namely NSs for which emission is believed to be powered 
by their strong magnetic energy. It is probable that the 
bursts were emitted by LS I +61°303 itself. Otherwise, 
we would be witnessing the unlikely alignment, within 
a couple or arcmin, of a gamma-ray binary (a popula- 
tion of objects for which a handful members are known) 
with a magnetar-like burst-emitting object (for which wc 
know 20 sources). If the LS I +61° 303 origin is accepted, 
any model of its multi- wavelength emission should thus 
provide an explanation of such bursts. 

Under the common assumption of pulsars emitting 
their rotational energy via magnetic dipolar losses, the 
NS surface dipolar magnetic field can be es timated from 
the observe d period and period derivative (|Pacinil [l967t 
iGoldl fl969() . For the known magnetars, this usually 
ranges from ~ 5 x 10 13 to 2 x 10 15 G; recently, how- 
ever, two sources with a lower field, ^7x 10 12 G, were 
discovered, the emission of which is still believed to be 
owered by non-di polar compone nts of the magnet ic field 
Rea et al.1 l2lH0ai iTurolla et al.1 [20TT1 : iRea et alJ[20Tl . 



About 20 magnetars are known to date, all being iso- 
lated pulsars with periods ranging from 0.3-12 s, usu- 
ally young spin-down ages ranging between 0.7-230 kyr 
(again with the two exceptions reported above which 
are also much older systems), and X-ray luminosities 
of the order of 1Q 33 ~ 35 erg s _1 (see iMereghettil 120081 : 
IRea fe Espositol2011l for recent reviews). Magnetars, his- 
torically divided into the two subclasses of Anomalous 
X-ray Pulsars (AXPs) and the Soft Gamma Repeaters 
(SGRs), display a large variety of bursts and flares, with 
properties at variance with those observed from other 
compact objects such as accreting NSs or BHs. Mag- 
netars bursts can be empirically divided in three main 
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TABLE 1 

Bursts observed by Swift-BAT from LS I +61°303 





Burst No.lM 


Burst No.Il( b ) 


Date 


2008 Sep 10 


2012 Feb 5 


Position uncertainty 


2.1' 


3' 


Ang. sep. 


0.60' 


1.07' 


Tioo (s) 


0.31 


0.044 


Fluence (10 _s erg cm -2 ) 


1.4 ±0.6 


0.58 ±0.14 


r 


2.0 ±0.3 


3.9 ±0.4 


Luminosity (10 37 erg s _1 ) 


2.1 


6.3 



Note. — The positional uncertainty is given at a 90% 
confidence level, including also systematic uncertainties. The 
angular separation is calculated with respect to the position of 
the optical counterpart. The Tioo duration and the fluences 
are estimated in the 15-50 keV band. Burst spectra were 
fitted by a power law with index T. The av erage luminosity is 
estim ated by assuming a distance of 2 kpc {Frail & Hi ellminel 

nam . 

M ITorres et al.l (2012) 

( b > From IBurrows et al.l J2012r ), see also 
http: / /gen. gsfc.nasa.gov/notices_s/513505/BA/ 

classes (although there is probably a continuum among 
them): the short bursts (~ 0.01 - Is; 10 37 " 40 erg s" 1 ), 
the intermediate bursts (~ 5 — 50 s; 10 40-42 erg s _1 ) and 
t he giant flares (~ 10 - 500s; 10 43 " 47 erg s" 1 ). 

ITorres et al.l (|2012t ) have started to study how a high- 
field NS could cope with the multi-wavelength phe- 
nomenology of LS I +61°303. In their scenario the NS 
would behave as a usual rotation-powered pulsar only 
when far from the companion star, whereas close to pe- 
riastron, the increased pressure exerted by the matter of 
the Be equatorial disk would rather overcome the pulsar 
pressure, quenching the rotation-powered pulsar behav- 
ior. Though, accretion of the matter captured would be 
inhibited by the qu ick rotation the NS, which w ould then 
act as a propeller (jlllarionov fc Sunvaevl[l975h . Such al- 
ternation between ejector and propeller states along the 
orbit, which we r efer to as a flip-flop st a te, w as origi- 
nally proposed by iGnusareva fc Lipunovl (|1985[ ) for NS 
in close-binary orbits of high eccentrici ty, and already 
applie d to the case o f LS I +61°303 bv iCampana et all 
(fl995h andE amano 

In this paper we delve further into this scenario, es- 
timating the interval of spin periods at which a NS in 
LS I ±61°303 is expected to behave either as an ejector 
or a propeller, and the duration of each of the states expe- 
rienced by a NS in an eccentric binary system, as it spins 
down during its initial evolution. By taking into account 
the constraints set on the parameters of the system by 
the observed 7-ray luminosity, we also estimate the rel- 
ative likelihood of observing the assumed NS in one of 
the different states. Results are discussed, comparing the 
case of an assumed high-field NS in LS I +61°303 to pos- 
sibly related systems, such as rotation powered sources 
in eccentric binary systems, as well as very long period 
HMXBs, thought to have host a magnetar in their early 
evolutionary stages. 

2. SPIN EVOLUTION OF A NS 

A NS evolves through different emission mechanisms 
during its existence, ejector, propeller, accretor, and 
georotator, depending on the balance between the out- 
ward pressure exerted by its electromagnetic field, and 



the ram pressure of the surrou nding matter (see, e.g., 
iLipunov et alJl992T : iGh osh 2007, and references therein). 
The electromagnetic pressure critically depends on the 
spin period of the NS, P, and on the strength of its dipo- 
lar magnetic field, B\. On the other hand, if the NS has 
a high mass companion, the pressure exerted by the mass 
lost by the latter, through a wind and possibly an equa- 
torial disc such as in the case of a Be star, is mainly 
determined by the density and velocity of the outflow, 
and by the velocity of the NS motion along the orbit. It 
turns out that, once the NS magnetic field and the rate 
of mass captured by the NS, Mi, are set, the state in 
which the NS lies is determined by its spin period, P. 
While at fast spin rates the NS behaves as an ejector, it 
is expected to become a propeller first, and subsequently 
accrete matter on its surface as it slows down. 

If the NS orbit is highly eccentric, the rate of mass cap- 
tured by the NS may vary by orders of magnitude along 
an orbital cycle, even if the companion star is assumed to 
lose mass at a constant rate; along its orbit the NS may 
then switch from one state to the other, such as it is the 
case for the flip-flop, ejector/propeller state proposed for 
LS I +61°303. 

2.1. Ejector state 

In the ejector state, a NS spinning at an angular fre- 
quency, il = 2-k/P, emits energy across the entire elec- 
tromag netic spectrum at the expenses of its rotational 
energy. ISpitkovskvl (|2006D evaluated the spin-down lu- 
minosity of a strongly magnetized oblique rotator, by 
solving for the dynamics of the field in presence of con- 
ducting plasma (the so-called force-free limit of rclativis- 
tic magneto-hydrodynamics): 

Q4 

L cj = (Bii? 3 ) 2 — (l + sin 2 a). (1) 

Here B\ is the dipolar magnetic field at the equator of 
the NS, Ri is the NS radius and a is the angle between 
the magnetic and the spin axis. The spin-down torque 
acting on the NS can be therefore expressed by 

iV cj = = -(5ii??) 2 ^J(l + sin 2 a). (2) 

According to t h e co nventio nal pulsar models 
(jGoldreich fc Julian! 119691: see ILipunov et ail 119921 
for a review), a NS behaves as an ejector as long 
as it manages to keep the surrounding plasma from 
penetrating into its light cylinder, the radius of which is 
Rlc = c/fi. To stop the in-fall of the matter captured 
by the gravitation of the NS before it penetrates into 
the light cylinder, the pressure exerted by the NS 
electromagnetic field must o vercome the pre s sure o f the 
in-falling matter. Following iBondi fc Hov lc (|1944f> . the 
radius at which matter is captured by the gravitational 
field of the NS is: 



where M\ is the NS mass and v re i is the velocity of the 
captured matter with respect to the NS. At lower radii, 
matter would start falling towards the NS at a velocity 
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of the order of the free-fall value, 



v ff 



2GMx 



exerting a pressure 



Pr 



2 _ Ml y / 2GM 1 
PV ff _ 4^ r 5/2 



(4) 



(5) 



where p is the gas density, and the mass continuity equa- 
tion was used. The pressure of the NS electromagnetic 
field outside the light cylinder, 



Alter 2 



(6) 



scales less steeply with the distance than the pressure 
of the incoming matter inside the gravitational radius 
(pram oc r~ 5 / 2 ); when the pressure of the captured matter 
evaluated at Rq overcomes the electromagnetic pressure, 
it is then expected to penetrate into the light cylinder as 
it falls inwards, driving the NS out of the ejector phase 
as a consequence. The matter in-fall may be stopped in 
fact only by the NS magnctosphcric pressure, 



Pu 



3\2 



(BiRt) 



1 



87rr 6 ' 



(7) 



In this case, the size of the magnetosphere is indeed de- 
fined in terms of the balance between p mag n and p iam , 
yielding the so-called Alfven (or magnetic) radius, 



Rm 



(BiRi) 



3\4/7 



M 2 / 7 (2GM 1 ) 1 /7 



(8) 



When the magnetosphere is able to extend up to the 
light-cylinder radius again (e.g., because of a decrease of 
the mass capture rate) , the NS may then resume to emit 
as an ejector. The condition to recover an unscathed light 
cylinder [R M > Rlc, i.e., Pram(-ffLc) < 
slightly different than the condition to stop the ejector 
mechanism [Pram(-RG) > Pc}(Rg)]- The former condition 
is fulfilled at a larger NS spin period when the other rel- 
ev ant magnitudes are held fixed (see, e.g. the discussion 
in iTorres et al.l |2012|) , and can be considered more re- 
strictive to identify when the NS abandons the ejector 
state. For simplicity and to be conservative, we then 
consider throughout this paper the equality, Rm = Rlc, 
to define the transition either from, and into the ejector 
state. The period at which the transition from the ejec- 
tor to the propeller state takes place can be expressed in 
terms of the rate of mass captured by the NS, once the 
the mass, radius and magnetic field of the NS are fixed: 



iM 1 ™-V7 12/7 



-2/7 



(9) 



sup prop (mi) = 0.24 b{ 

h = (Bi/10 13 G), mi = (Mi/1.4 M ), n = 
(Mi/10 17 g s _1 ) are the magnetic 



Here 

(i?i/10km) and rhi 
field, the mass, the radius, and the rate of mass captured 
by the NS, in units of the values we consider as fiducial 
in the rest of the paper (see Table [2] for a complete list 
of the scale units considered) . 

2.2. Supersonic propeller state 

When the NS stops acting as an ejector, accretion onto 
its surface is still inhibited by the rotation of the NS. The 



TABLE 2 

Scale units used in this paper. 



Scale 



Definition 



hi 


Bi/10 13 G 


NS magnetic field 


'"i 


Mi/1.4 M 


NS mass 


ri 


Rl/lOkm 


NS radius 


1 


il/10 45 g cm 2 


NS moment of inertia 


iii i 


Mi/10 17 g s" 1 


NS mass capture rate 


m™ ax 


Mr°7l0« g s- 1 


max NS mass capture rate 


m™" 


Mr m /5 x 10 12 g s- 1 


min NS mass capture rate 


b 2 


B2/O.6 kG 


Be star magnetic field 


in-2 


M 2 /12.5M Q 


Be star mass 


r> 


K2/IO Rq 


Be star radius 


H-2 


n/2 


index of Be disc mass capture 




rate radial dependence 


<h 


d C ut/7R2 


Be disc cut-off size 


■ P 

m\ 


Mf/10 18 g s" 1 


Be star mass loss rate 


V 


■"So/10 8 cm s- 1 


Be star wind term, velocity 



magnetosphere spins much faster than the in-falling mat- 
ter at the boundary defined by Rm, and the interchange 
instabilities allowing the the plasma t o enter into the 
magn etosphere are strongly suppressed (|Elsner fc Lamb! 
[1973). To express the centrifugal inhibition of accretion, 
we introduce the co-rotation radius, defined as the ra- 
dius at which the linear velocity of the rotating magneto- 
sphere equals the Keplerian rate, H^(r) = (GA/i/r 3 ) 1 / 2 : 



R c 



( GM\ 



\ ft 2 



1/3 



(10) 



and define the NS fastness as ([Ghosh fc Lamblll979| ). 



» = (Rm\ 

n K (R M ) \RcoJ 



3/2 



(11) 



If Rm > Rm (i.e. uj* > 1), a centrifugal barrier prevents 
the accretion of matter onto the N S surface, and the NS 
is said to lie in a propeller state. Illlarionov fc Sunvaevl 
(|1975l) expressed the luminosity emitted by the NS in 
this state in terms of the energy needed to balance grav- 
itational energy of in-falling matter, £pf op = Mf 2 f / 2. 



Subs e quent studies (iDavies et al.lll979l ; lDavies fc Pringlej 
[l98lt iMineshige et alJ 11991 " showed how a quasi-static 
corona forms around the NS, as far as the quickly ro- 
tating magnetosphere transfers energy to the incoming 
matter at a rate larger than the gas cooling rate. Such 
a corona extends down to the radius where its pres- 
sure is balanced by the pressure of the magnetic ficlcQ, 
Rin ~ Rm- At the interface between the corona and 
the magnetosphere the gas is shocked by the supersonic 
motion of the field lines, and energy is transferred to the 
coronal gas through turb ulent or convective motions (see 
iWang fc Robertsonlll985l for a detailed treatment). Such 
a transfer takes place at the expenses o f the spin of the 
NS, which decelerates at a rate (see, e.g. . IMineshige etaD 

1 A proper assessment of the magnctosphcric boundary depends 
on the details of the structure of the corona. However, the ra- 
tio between Ri n and Rm evaluated by a number of authors (e.g. 
IDavies fc Pr irmlc 1981 who evaluated it as (Rg/Rm) 2/9 ) is of the 
order of one for the parameters considered here. In fight of the 
large uncertainties on the plasma capture process (see we then 
consider Ri„ = Rm to make simpler the evaluation of the torques. 



4 



A. Papitto, D. F. Torres, & N. Rea 



[Mil): 



4axR 2 M v t (R M ), 



(12) 



where e is the energy density transferred by the NS dur- 
ing each revolution, vt is the velocity of the turbulence 
developing at Rmi and inRj^Vt (Rm) is volume of the 
gas perturbed by the motion of the dipolar magnetic 
field, assumed to be tilted with respect to the spin axis. 
As long as the linear velocity of the magnetosphere ex- 
ceeds the sound speed at Rm (taken to be of the or- 
der of the free-fall velocity, «//), and the energy re- 
leased to the corona by the NS domin ates radiative losses 
(see the discussion ofllkhsanov 2002, relative to wind-fed 
close binary systems), the propeller is considered super- 
sonic and the turbulent motions take place at a velocity 
Vt — Vff. Though, there is no gene ral consensus on the 
esti mate of the prope l ler effi ciency. iDavies et al.l (|1979f ) 
and lDavies fc Pringld (|1981f) consider e ~ pv 2 f f/ 2, recov- 
ering the scaling of llllarionov fc Sunvaevi ([1975D 



N, 



is 



M ^* 



(13) 



On the other hand, if e ~ pjilsR iyt) 2 /2 is considered, a 
much strong er torque is obtained (jMineshige et al.lfl991t 
IGhoshl[T99a . 



N, 



prop 



-\m^Jgm 1 r m 

6 



(14) 



where the numerical factor takes into account the de- 
gree of non-axisymmetry of th e magnetosphere with re- 
spect to the spin axis (see also [Wang fe Robertsonl | 1985 |: 
llllarionov fe Kompaneetsl IT9 90; Bis novatvi-KoganTll99lL 
who derived expressions with a similar scaling) . The dif- 
ference between the two estimates given by Eq. (fl~3j) and 
([14)) may be as large as ~ 10 4 , when a NS with the fidu- 
cial parameters defined in Table [2] is rotating close to the 
critical period marking the transition between the ejector 
and the propeller state (Eq. [9]), since cj* « 100 in that 
case. The discrepancy between the two propeller torques 
yields a significant uncertainty in the evaluation of the 
timescale of the NS ev olution in the propeller state (see, 
e.g., the discussion of iFrancischelli fc Wiiersl 120021 and 
IMori fc RudermanI I2003D . In order to be conservative 
when describing the evolution of LS I +61°303 in the 
flip-flop state, we consider the two estimates given above 
as limiting cases. 

2.3. Subsonic propeller and onset of accretion 

As the velocity of the NS decreases and becomes com- 
parable to the speed of sound at the inner boundary 
of the corona, the propeller becomes subsonic, with the 
turbulence t raveling at a speed set by the NS rotation, 
v t ~ QRm ([Davies fc Pringld 1198 1). We assume that 
such a transition takes place when Rm — Rco, translat- 
ing into a period: 



sup prop— ^sub prop 



18 b\ 



6/7 . -3/7 -5/7 18/7 



(15) 



From Eq. (|12p it is deduced how the torque experienced 
by the NS in this stage differs by a factor CIRm /v / /(Rm) 
from those defined in the previous section. The rate 



at which energy is transferred from the NS to the sur- 
rounding corona decreases with increasing period; gas 
in the corona then starts to cool down, facilitating the 
plasma entry into the magnetosphere. A fraction of the 
incoming matter may then accrete down to the NS sur- 
face already in the subsoni c propeller state (the so-called 
settling regime studied bv lShakura et al.ll2012l ). Subse- 
quently, the cessation of any significant barrier effect is 
achieved when the energy released by the rotating mag- 
netosphere to the incoming matter can be neglected with 
respect to the cooling of the gas; the spin period for such 
a transition was estimated as: 

-Psub prop^acc = 91 bf 121 m^ 5/T m^ 4/21 r{ 6/7 s, (16) 

bv llkhsanovl (|2001D . who used this expression to estimate 
the duration of the subs onic propeller s tate in wind-fed 
binary systems (see also llkhsanov|[2007| ) . 

3. MAGNETIC DISSIPATIVE TORQUES 

The interaction between the strong magnetic field of a 
magnetar and the field of a low mass, convecti ve compan- 
ion star was invoked bv lPizzolato et all ([2008) to suggest 
how the spin period of the X-ray source IE 161348-5055, 
could have been locked to the orbital period of the sys- 
tem, similar t o what happe ns in polar cataclysmic vari- 
ables (see e.g. IWarner| [l995;). In such a case, it is in fact 
argued that a dissipativc torque, 



d 3 ' 



(17) 



develops as the magnetic dipole of the white dwarf, fj,i, 
and the magnetic field of the companion star, p 2 , inter- 
act at an orbital separation d. The magnetic field of the 
companio n star can be either induced by the white dwarf 
field fe.g. [Toss et al.l H971: lLamb et aL|[T983t ICampbelll 
19841) or intrinsic to t he companion star ( Campbeillll985l : 
Hameurv et al.|[l987t ) . The rotation of a low mass com- 
panion star belonging to a close system (P or b ~ few 
hours) is synchronized to the orbital motion by tidal 
forces on a relatively short timescale, « 10 2 — 10 3 yr, as 
it is o btained by considering the relation given bv lZahnl 
(1977) 



1 (M 2 



6 V Mi 



M 2 Rl 



1/3 



k 2 M 2 R\ \R 2 



(18) 



Here, M 2 , R 2 , L 2 and I 2 are the mass, radius, luminosity 
and momentum of inertia of the non degenerate star, k 2 
is the consta nt of apsida l motion and is of the same order 
of I 2 /M 2 R\ (IZahnlll977D . and a is the semi-major axis of 
the orbit. The torque due to the interaction between the 
magnetic fields of the two stars then acts to bring the 
white dwarf to synchronicity with the orbit, as well. 

The projected rotational velocity of th e ~ 12.5 M ff 
Be sta r in LS I +61°303 was measured bv lCasares et al.l 
(|2005| ) as 113 km s _1 . Such a value corresponds to a 
spin period of ~ 4.5 r 2 smi d, where i is the inclina- 
tion of the system and r 2 = (R 2 /\QRq) is the radius of 
the companion in units of 10 R@. Plainly put, the ro- 
tation of the Be star is not locked to the orbital period 
of the system. This is consistent with the time needed 
to synchronize the spin of the non-degenerate star to the 
orbit of LS I +61°303 through tidal interactions with the 
compact object (w 10 7 yr , see Eq. (fl"8[0 . 
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However, the magnetic field of a star with the proper- 
ties of the Be star in LS I +61°303 can be in principle 
very intense. A surface field of ^0.6 kG wa s measured 
from a star of the same luminosity class (jPetit et al.1 
US!), while even larger fields were measured from pe- 
culiar sta rs with spin periods lower than 1 d (see e.g. 
table 1 in lOskinova et all 1201 lL and references therein). 
An intrinsic field of this order is larger by orders of mag- 
nitude than any field which may be induced by the NS 
field, and the magnetic dipole moment of the Be star, 
ps 10 38 G cm -3 , would be much larger than the NS mo- 
ment, ps 10 31 G cm~ 3 , as well. However, in the case of 
a relatively wide binary such as LS I +61° 303, the steep 
dependence of Eq. (jlTf on the orbital separation greatly 
reduces the magnitude of the torque. 

4. MASS CAPTURE 
The rate of mass lost by a Be star such as that in 
LS I +61°303, Mi, is given by the sum of the mass lost 
through a fast polar wind, Mf , and a slow equatorial disc, 

M% (e.g. iWaters et all 119881 ). According to the Bondi- 
Hoyle description, the mass lost by the companion star is 
then captured by the NS at a rate, M x = M 2 (R G / d) 2 / 4, 
where Rq is the radius of gravitational capture defined 
by Eq. Q. The estimate of M\ obtained under this 
approximation is thus very steeply dependent on the ve- 
locity of the captured mass with respect to the NS, V re i. 
The velocity of the polar, radiatively driven outflow is 
described by v p (r) ~ ^(1 — R^/ d), where v 1 ^ is the 
polar wind terminal velocity (e.g. lLamers fe Cassinellil 
Il999f ). The orbital velocity of the compact object in 
LS I +61°303 can be n eglected with respe ct to the wind 
velocity (see Fig. 13 of lTorres et al.|[2012h . and the rate 
of mass captured from the polar wind is given by: 



Mf 



2 (vL)W 



Ih 
d 



(19) 



respec- 



In the following we scale the rate of mass lost by the Be 
star and the wind terminal velocity in units of m p 

(Mf/10 18 ) g s- 1 and v = «,/10 8 cm s" 
tively, of the order of the estimates given bv lWaters et al.l 
(jl988fl . 

A rapidly rotating Be star loses mass through an 
equatorial disc at a rate which is generally 10 — 100 
times larger than that of the polar wind (see, e.g., 
lLamers fe Waterslll987D. By m odeling the observed IR 
excess. iMarti fe Paredesl (|1995[ ) estimated the equatorial 
disk mass loss rate of the Be star in LS I +61°303, M^ , to 
lie between 0.25 and 2.5xl0 19 g/s. This range depends 
on the value of the radial velocity of the flow at the sur- 
face of the Be star, assumed to vary in the range 2-20 
km s _1 . However, modeling the capture of mass from 
the equatorial disc is much more uncertain with respect 
to the polar wind case. In this case, for any reasonable 
assumption on the radial and azimuthal velocity profile 
of the equatorial disk matter, the orbital velocity of the 
compact object in LS I +61°303 cannot be ignored, as it 
is the case when the poloidal wind is considered. Most 
importantly, the Bondi-Hoylc approximation yields val- 
ues of the mass capture rate at the periastron of the or- 
bit which may be unphysical larger than the rate of mass 
lost by the Be star, essentially because of the low rela- 



tive velocity of the disc matter in such a system ( ;$ 10 7 
cm s _1 ; see Fig. 13 in lTorres et al.ll2012[ ). Given these 
large uncertainties, to model the capture of mass from 
the equatorial disc we parametrize with a power law its 
dependence on the orbital separation, d, adding a strong 
cut-off at a distance d cut to reproduce a tidal truncation 
of the disc due to the interaction with the compact object 
(jOkazaki et al.ll2002t ). We thus consider: 



Mf = 



exp 



d 

dcut 



10 



(20) 



Here, M™ ax = Mf{d m i n ) is the maximum rate of mass 
captured by the NS from the equatorial disc, when the 
NS is close to the periastron of the orbit, i.e., at an 
orbital separation d m i n = a(l — e). In the following, 
we consider the eccentr icity e = 0.63 ± 0.11 measured 
by iCasares et al.l (I2005D. and comp atibl e with the es- 
timate s of iGrundstrom et all (|2007f ) and lAragona et al.l 



([2009h . while the semi-major axis of the orbit, a, is es- 
timated as 6.3 x 10 12 cm by using the third Kepler law 
for a system with an orbital period of 26.5 d and a to- 
tal mass of 14 M Q . Since the disc is much denser than 
the polar wind, the maximum rate of mass captured by 
the NS, M™ ax is equal the rate of mass captured from 
the disc, M 1 ' max ; in the following we scale this value 

in units of rh" ax = M i maa /10 17 g s _1 , which was also 
used by iDubud (|2006f ). Such a value is roughly in be- 
tween the estimates considered by IZamanovl ( 19951 ) and 
IGreeorv fe Neishl ((20021) . ps 3 x 10 17 and ps 0.6 x 10 17 g 
s _1 , respectively. A valu e of the same o r der, ~ 0.5 x 10 17 
g s _1 , was also found bv lRomero et a l. (2007) from sim- 
ulations of the interaction between the equatorial disk 
of the Be star in LS I +61°303 and the compact object, 
assumed to be accreting in the case they considered. It 
is also worth to note how a peak accretion rates up to 
~ 6 x 10 17 g s -1 , is deduced from the X-ray luminos- 
ity observed at the peak of the outburst sho wn by the 
pulsar, 4U 0115+63 (see iFerrigno et al.|[20TlT and refer- 
ences therein), which has a B0.2Ve star companion and 
orbital parameters (P or b = 24.3 d, e = 0.34) similar to 
those of LS I +61°303. The index of the power law n 
of Eq. (|2U|) is varied between 1 and 3, to qualitatively 
reproduce the dependence of the mass c apture rate on 
the di stance found by the simulations of iRomero et al.l 
(|2007f ). and in particular the ratio in the range 10-100 
between the maximum and minimum mass capture rate 
in the absence of a significant cut-off (see the solid curve s 
in the left panel of Fig. 13 plotted by iTorres et aLll2012t) . 
The fiducial unit of the truncation radius of the equato- 
ria l disc is set to 7 R 2 , of the order of the estimates given 
bv lGrundstrom et~a l. ( 2007) on the basis of the observed 
equivalent width of the Ha emission li ne. Larger values 
are a lso reported in literature (see e.g.. lGregorv fc Neishl 

HQQl who give ~ 12 R 2 ). 

To evaluate the orbital dependence of the wind and disc 
contributions to the total rate of mass captured by the 
compact object, Mi = M^ + Mf , we use the relation d = 
a(l — e cos e) to express the orbital separation in terms of 
the eccentric anomaly, e. A red solid line shows in Fig. [T] 
the mass capture rate as a function of e, for the fiducial 
set of parameters. In the following we shall also consider 
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0.4 0.6 
Eccentric anomaly 

Fig. 1. — Rate of mass captured by a NS from the equatorial 
disc (blue dotted line) and the poloidal wind (magenta dotted line) 
emitted by the Be star in LS I +61°303, as a function of the eccen- 
tric anomaly, and for the fiducial values of the relevant parameters 
(sec Table [2J . Red solid line is the sum of these two contribu- 
tions. The black dashed line shows the case of an increased Be star 
mass loss rate, rh'[ iax = 5, and a disc cut-off beyond the maximum 
orbital separation, d cu t > a(l + e). 

a mass capture rate increased by a factor of five, with an 
outer radius exceeding the maximum orbital separation 
(see black dashed line in Fig. []}, in order to mimic an 
enhanceme nt of the rate of mas s los s of the order of that 
reported bv lGregorv et al.|[l989l and lZamanov et al.lll999l 
to explain the observed super-orbital variability (see e.g., 
lGregory||2002l ). with the possible expulsion o f mass shells 
in the disc (see, e.g. lGregorv fc Neishl 120021 ). 

5. TIMESCALES AND CHANGES OF STATE 

The time it takes for a NS to reach a period P, under 
the action of a spin-down torque, N(P), is obtained from 
the integration of the equation, 

2-kI dP 

between the initial period, Po, and P: 

dP 



N(P) = m 



(21) 



t 



-27r/ 



Po P 2 N(PY 



(22) 



A NS is generally considered to spin at its birth at a 
period of few tens of ms, therefore spinning down as a 
rotation-powered pulsar. If it belongs to a binary system 
such as LS I +61°303, the ejector phase will end as the 
pressure of the mass captured at periastron overcomes 
the pulsar pressure. The system then enters in the flip- 
flop state (i.e., the state at which it is in ejector phase in 
apastron and in supersonic propeller in periastron) when 
its period attains a value -P e j->flipftop- This is obtained by 
evaluating Eq. ^ (plotted as a blue solid line in Fig. [21 
where the periods at which the various state transitions 
take place are plotted as a function of the mass capture 
rate, for a magnetic field b\ = 1), at the maximum ac- 
cretion rate experienced by the NS along its orbit, m r [ MX 
(rightmost vertical dashed line in Fig. 0): 



P 



cj— ^fiipflop 



n o/i u 4 / 7 -V 7 12 / 7 
0.24 b x m 1 r x 



f-max\ 

(m 1 ) 



-2/7 , 



(23) 
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Fig. 2. — States accessible to a system like LS I +61°303 in the 
NS spin period vs. mass capture rate phase space. Solid lines mark 
the transitions between the ejector, supersonic propeller, subsonic 
prop eller, and accretor states, as defined in text [see equations l|9)l . 
115> . and 116H . and evaluated for 6i = 1. Vertical dashed lines 
mark the minimum and maximum capture rate experienced by the 
NS along its orbit, while the horizontal dashed lines indicate the 
values of the spin period at which transitions take place, when 
the relevant parameters are set equal to their fiducial values (see 
Tabic EJ). 



Assuming for the moment the magnetic field of the NS 
to be constant, the duration of the ejector phase is found 
by evaluating Eq. (|22|) between Po and P j->fiipflopj with 
the ejector torque being described by Eq. ([2]): 



t ej =6.3b- 6/7 (m r l iax ) 



- 4 / 7 ™- 2 / 7 T.- 18 / 7 



1(1 



sin 2 a) 



_1 kyr, 
(24) 



where I = (Ji/10 45 ) g cm 2 and I\ is the NS moment of 
inertia. The value of t e j does not significantly depend on 
Po as long as it is much smaller than P j->.ni p fiop- 

Evaluated at a given period, Eq. ([9]) implicitly defines 
a critical value for the mass capture rate at which the 
NS switches from the ejector state to the propeller state 
(and vice versa), M cr it(P). By using the relations given 
in the previous section to express M\ as a function of 
the orbital separation, this condition translates into a 
critical separation d cr u(P) 1 defining the portion of the 
orbit during which the NS behaves as a propeller [for 
d < d cr it(P)], or as an ejector [for d > d cr i t (P)]- We 
then express the torque experienced by the NS during 
the flip-flop state as: 



N s (d) = { ^ 

lyprop \ iy magn 



for d > d crit (P) 
for d < d crit (P) 



We summarize in Table [3] the spin evolution timescales, 
t = Q/Q, implied by the torques introduced in the pre- 
vious sections. The average torque experienced by the 
NS along one orbit is: 



< N ff > 



1 



Porb 

l 

Pnrh 



N s (t)dt 



2tt 



N s [d(M)} 



dt 



dM 



dM. (25) 
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TABLE 3 

Spin evolution timescales 



Torque 



Timescale, Q/Q (kyr) 



EM 

IS Supersonic Propeller 
G Supersonic Propeller 
IS Subsonic Propeller 
G Subsonic Propeller 
Magnetic Torque 



144 &i r~ b XP 2 
77.54 b 4 / 7 m~ 9/7 m~ 



1.43 b~° 

a i I.-2/7 
0.1 0, ' m 

-2 



-3/7 , 
"1, 

-6/7, 



8/7 12/7 T 
' r 1 IP 

2/7 



-24/7 . 



3/7-6/7 



1 "'1 ' 1 
0.11 b~ z mi r~ b IP 

21 x 10 3 6i rf f> 2 rf d~ 3 2 P" 1 



1P- 



Note. — See Table [2] for the definition of the scale units. A 
value of a = 45° was considered to estimate the electromagnetic 
torque. P is the NS spin period in seconds. 

Here, Ai is mean anomaly, \dt/dM\ = P or i,/2ir, the 
relation between the mean and the eccentric anomaly is 
given by the Kepler's equation, M = e — e sin e, while 
the orbital separation is obtained as d = a(l — ecose). 
The magnitude and dependence on the orbital separation 
of the different torques introduced so far, as well as the 
critical distance at which the transition ejector/propeller 
takes place, evaluated for example at a period P = 1 s 
and for the fiducial parameters defined in Table [2j are 
plotted in Fig. [3] 

The NS abandons the flip-flop state when it reaches 
a period, -Pfiipflop^supprop, such that it is in a propeller 
state at all points in the orbit, even when the rate of 
mass capture is minimum, i.e., at apastron, M cr n ~ 
M\{d max ). It can be seen from Fig. [T]that, as long as 
the equatorial disc is cut off at a distance d cu t < d max , 
this value is set by the rate of mass captured by the 
polar wind, M^{d max ). Considering a scale unit of 
m™" = M™ m /5 x 10 12 g s _1 for the minimum rate 
of mass capture, of the order of that expected at the 
apastron of the orbit if the Be star loses mass through 



the polar wind at a rate of 10 18 

-Pflipflop 



Pi 



fiipflop— ^sup prop 



p prop is found from Eq. ([9]) : 

= 4.1 b\ n (rh? in )- 2 / 7 : 



g s 1 , the value of 



-1/7 12/7 

1 ' V ' 



(26) 

Fig. [2] shows how for the typical values of the parameters 
relevant to the case of LS I +61°303, the NS leaves the 
flip-flop state well before its propeller torque becomes 
subsonic, and only the torques defined in § !2.1l and § 12.21 
are relevant to the evaluation of the time spent by the 
NS in the flip-flop state. 

6. RESULTS AND CONSTRAINTS 

6.1. The flip-flop phase 

The total time spent in the flip-flop phase by a NS 
in LS I +61°303, is evaluated by integrating Eq. ([22]) 
between Pej->fiipflo P and Pfii P fio P -«upprop- For the set of 
fiducial values we considered, the source stays in a flip- 
flop state when its spin period lies in the range between 
0.24 and 4.1 s. The total time spent in this state cru- 
cially depends on the relation considered to express the 
propeller torque. Values of 25 and 282 kyr are obtained 
when the torque of Eq. ([M]) and (fT"3|) are considered, 
respectively. Such a large discrepancy is due to the dif- 
ferent dependence on the system fastness of the two 
expressions. On the other hand, the torque due to the 
interaction between the magnetic fields of the NS and 
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Fig. 3. — Propeller torque evaluated according to the two cases 
pres ente d in text (red and magenta solid curves, respectively, see 
Eq. I13I and I14I ). eject or torque (blue; see Eq. [2]) and magnetic 
torque (cyan; see Eq. I17M . plotted as a function of the orbital 
separation between the two stars of LS I +61° 303 (expressed in 
units of the semi-major axis), evaluated for values of the relevant 
parameters set equal to the fiducial units (see Table [2j , and for 
an example period of P = 1 s. Vertical dashed line mark the 
critical separation, d cr it(P = Is), at which the system performs 
a transition from ejector (d > d cr it(P = Is)) to propeller (d < 

dcrit(P = Is)). 

of the companion star (see § [3]) has a negligible impact 
on the evaluation of the flip-flop time-scale even if the 
weakest propeller torque is considered. The timescales 
obtained have to be compared with the interval of rj 10 3 
kyr it would take for a NS with the assumed fiducial val- 
ues to cover the same interval of periods, only by spin- 
ning down as a rotation-powered pulsar with a constant 
magnetic field. 

In the following, we focus on the results obtained with 
the stronger propeller torque, Eq. (|14l) , since it gives val- 
ues which can be considered as conservative lower limits 
on the total time that a NS in a system like LS I +61°303 
is expected to spend in the flip-flop phase. We plot in 
Fig. @] the evolutionary track of the NS spin obtained by 
considering the fiducial parameters defined above. The 
dependence of the total time spent by the system in the 
flip-flop state can be approximated as 



-0.3 .-1.1 



V 1 « 



r l, ^kyr. (27) 



Here, d*? = {dcut/TRi) and n% = in/2). The flip-flop 
timescale depends strongly on the strength of the NS 
magnetic field, and on the amount of mass captured by 
the equatorial disc of the Be star, as it is expressed by 
the dependence on its size, d cu t, and less strongly on the 
maximum mass capture rate, m™ ax . 

Values of the magnetic field in excess of 2 x 10 14 G 
reduce the total flip-flop timescale to less than a kyr. 
Both, the range of periods for which the NS is in the flip- 
flop state and the magnitude of the spin down torque, 
increase when a stronger NS magnetic field is consid- 
ered, but the latter to a larger extent. The total flip- 
flop timescale also depends significantly on the amount 
of mass captured by the equatorial disk of the Be star. 
Such a time varies between «50 and 15 kyr when the 
maximum rate of mass captured varies between 10 16 and 
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Fig. 4. — Evolution of the spin period of a NS in a system such 
as LS I +61°303, obtained considering the fiducial values of the 
system parameter (see Table [2). Dashed horizontal lines mark the 
transition among different states, while the gray shaded area mark 
the time interval during which the NS is expected to lie in the flip- 
flop state. The black filled circle marks the spin period at which 
the system emits an ejector luminosity of 5 X 10 35 erg s _1 . 

5 x 10 17 g s" 1 ; a range of 17-36 kyr is obtained when 
the size of the disc cut-off takes value between 10 and 5 
i?2, while a smaller variation of « 15% is obtained when 
values of n in the range 1-3 are considered. 

6.2. B-field evolution 

We have also studied the effect of the decay of the dipo- 
lar magnetic-field on the total time spent by the system 
in the fl ip-flop state, by consi dering the simple relation 
given bv lAguilera et all (|2008l ) 

B l {t)=B Q — , e *?~ t/T ° , I; ft +S,., (28) 
1 + {toJth)\}- - cxp(-t/T )J 

where th and to are the timcscales for Hall and Ohm de- 
cay, set equal to 1 and 10 3 kyr, respectively, and B$ and 
.B as are the initial and asymptotic value of the magnetic 
field, respectively The evolutionary tracks evaluated for 
-B as = 5 x 10 12 G and a number of values of B are plotted 
in Fig. [5] For the larger initial field values we considered 
(50 x 10 13 G), the time it takes for the NS to enter the 
flip-flop state is so short (~ 0.2 kyr) with respect to the 
assumed value of th , that no significant field decay has 
still taken place and the timescale of the flip-flop state 
are correspondingly very short ( ,$ kyr). On the other 
hand, values in excess than 10 kyr are spent by the NS 
in the flip-flop state if the initial field is ;$ 5 x 10 13 G. 

6.3. Constraints from the apastron luminosity and 
Be-star mass loss rate variations 

We plot in Fig. [5] the different states (ejector, flip-flop, 
propeller) in which a NS in LS I +61°303 is expected 
to lie during the early stage if its evolution, depending 
on the values of its magnetic field and spin period; the 
tracks delimiting the various regions were evaluated from 
Eq. ()23[) and (|26[) , for fiducial units of the maximum and 
minimum mass capture rate. 

A relevant constraint on the possible position of 
LS I +61°303 in such magnetic field vs. period diagram 




0.1 1 10 100 

Time (kyr) 



Fig. 5. — Evolution of the spin period as a function of the time 
elapsed since the NS birth with Po = 0.01 s, for values of the 
initial NS magnetic field equal to (from top to bottom) 50, 10 and 
5 X 10 13 G. From left to right along each track, black points refer 
to the ingress in the the flip-flop state, the period at which the 
system emits an ejector luminosity of 5 X 10 erg s — 1 , and the 
egress from the flip-flop state. Numbers along each track denote 
the strength of the magnetic field at the ingress and egress from 
the flip-flop state, respectively, in units of 10 13 G. 

can be drawn on the basis of the flux observed when 
the NS is close to apastron. According to the simplest 
flip-flop scenario discussed (i.e., not considering the pos- 
sible effect of a relativ istic wind on the sp in-down of a 
magnctar, discussed bv lHarding et al.lll999t) . when in the 
apastron region the NS is powered by the ejector mech- 
anism, and its luminosity cannot exceed the value given 
by Eq. @. By evaluating such a relation for a = 45°, 
and an ejector luminosity of 10 37 , 10 36 , and 10 35 erg s" 1 , 
the tracks plotted as red solid lines in Fig.|6]arc obtained. 

Excluding the contribution of the bright Be star, the 
broad-band spectral energy distribution of LS I +61°303 
peaks in the MeV-Ge V band (se e, e.g., Fig. 6 in 
Chernvakova et ai]|2006l and Fig. 2 in lGupta k, Bottcherl 
2006J, and references therein). The 0.1-300 GeV flux ob- 
served by Fermi Large Area Telesc ope when the NS is at 
apastron was recently estimated bv lHadasch et al.1 (|2012l ) 
as 5.1(3) x 10~ 10 erg cm -2 s _1 , corresponding to a lumi- 
nosity of ~ 2.5 x 10 35 d 2 . erg s _1 , where 6.2 is the distance 
to the source in units of 2 kpc. Such a value is ~ 20% 
lower than the flux observed when the compact object is 
close to periastron. A flux of ~ 2 x 10 -10 was observed 
by COMPTEL on-board the Compt on Gamma-Ray Ob - 
servatory, in the 1-10 MeV band (jTavani et al.l I1996Q . 
while the source is significa ntly dimmer in the X -ray and 
TeV energy bands (see, e.g. lHadasch et al"1l2012l and ref- 
erences therein). 

If such luminosities are ejector-only generated, the 
spin-down power of the NS must lie in a range between 
few x 10 35 and few x 10 37 er g s _1 (see, e.g., the dis- 
cussion of Z abalza et al.l 120 111 who considered an ef- 
ficiency of spin-down to 7-ray luminosity conversion of 
0.034(L ej /10 36 erg s^ 1 )" 1 / 2 , following lAbdo et al.ll2uToft : 
the lower end of this interval is obtained when the emis- 
sion is beamed in a 1 sr solid angle, while the higher val- 
ues correspond to more isotropic pulsar beaming models. 
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Fig. 6. — Ejector, flip-flop and propeller states plotted in the NS 
magnetic field vs. spin period phase space, evaluated for a NS in 
LS I +61° 303 and for the fiducial values adopted for the maximum 
and minimum mass capture rate (m max = 1, fa mln = 1). From 
top to bottom, the red solid lines mark the relation between the 
period and the magnetic field of the NS when the ejector luminosity 
is 10 37 , 10 36 and 10 35 erg s — , respectively, and the magnetic offset 
angle is a = 45° . 

It is then clear how a pulsar with such an ejector luminos- 
ity is unlikely to lie to the right of the red line labeled as 
10 35 in Fig. |6l as well as to the left of the line labeled as 
10 37 . Even recalling how the displacement of the tracks 
plotted in Fig. [5] depends on the exact values of unknown 
parameters, such as the maximum rate of mass captured 
by the NS, as well as on the assumed ejector braking law 
and details on the ejector/propeller transition, the rela- 
tively large gamma-ray luminosity indicates that the NS 
must be young, and close to the transition between the 
ejector phase and the flip-flop state. 

Thus, if the source lies in a flip-flop state, its period 
must be close to the value it had at the beginning of 
such a phase, i.e. that given by Eq. (f23|) . This obvi- 
ously reduces the time spent in the flip-flop state before 
reaching a period corresponding to a certain ejector lu- 
minosity (for a given magnetic field), with respect to that 
given by Eq. (|27|) . For the case of a maximum mass cap- 
ture rate of m\ MX = 1 and a magnetic field b\ = 1, the 
time spent in such a state before spinning down to a pe- 
riod corresponding to an ejector luminosity of 5 x 10 35 
erg s _1 , is 3.1 kyr. By considering larger values of the 
ejector luminosity, the range of periods for which the NS 
is in the flip-flop phase further reduces (see Fig. [7J|); a 
timescale of 1.8 kyr is obtained for an ejector luminosity 
of 10 36 erg s _1 , while extremely small values, ,$ 0.1 kyr, 
are obtained if L e j £ 2.5 x 10 36 erg s _1 . From Fig. |JJ] 
it is also clear that the interval of periods for which the 
ejector luminosity deduced from the observed 7-ray flux 
is compatible with a flip-flopping behavior increases for 
lower magnetic fields and larger maximum mass accre- 
tion rates, with respect to the fiducial values considered. 

Alternatively, the luminosity observed when the source 
is close to the apastron has bee n interpreted as powered 
by the propeller emission fe.g.. lBednarekll2009l ). In such 
a case LS I +61°303 would lie at the right of the flip- 
flop/propeller transition (green area of Fig. [6]). How- 
ever, even considering the stronger propeller torque in- 



troduced, Eq. (|14p . the expected luminosity is: 

r _ o Af G ~ 
-^prop — " iv prop — 

~ 0.9 x 10 36 b\' 7 m\ /7 rf' 7 m~ 2/7 P~ 2 erg s" 1 (29) 

where P is the NS spin period in seconds. Such a lumi- 
nosity is then compatible with the luminosity observed 
from the system when the NS is close to periastron (i.e. 
for mi ~ 1), if the period is £ 1 s, i.e. if the system 
is close to the transition between the ejector and the 
flip-flop phase. On the other hand, if the spin period is 
^ 4 s, as it is implied by the assumption that the source 
is always in the propeller state, such an effect is barely 
capable to power the observed emission close to perias- 
tron (riii ~ 1), for a NS magnetic field £ 10 14 G (i.e. 
61 }t 10). Close to apastron instead, the rate of mass 
capture is much lower, m ~ 10~ 4 , and the propeller falls 
short by ~ two order of magnitudes in accounting for the 
observed 7-ray luminosity (but see below). 

While the timescales derived so far do not depend 
much on the value of the minimum mass capture rate, 
m™" - , this parameter plays a major role in determining 
the state of the NS at orbital phases close to apastron, 
and therefore the value of the period at which the sys- 
tem abandons the flip-flop state (i.e. the location of the 
track labeled as flip-flop /propeller in Fig.[fj]). In particu- 
lar, an increase of the rate of mass captured by the NS, 
possibly caused by an increase of the rate of mass lost 
by the Be star and/or by a growth of the size equato- 
rial disc, may bring the system out of the flip-flop state. 
The magnitude of the variation of the mass capture rate 
needed to determine such a transition obviously depends 
on the initial location of the source in a diagram like that 
of Fig. UJJ If the source lies close the flip-flop /propeller 
transition, a variation by a factor of few of the mass cap- 
ture rate at apastron, owing for instance to an enhance- 
ment of the rate of mass lost by the Be star through its 
polar wind, is sufficient to inhibit the ejector emission 
throughout the orbit. On the other hand, if the source 
lies close to the ejector / flip- flop, as it is deduced from 
the observed 7-ray luminosity, a much larger increase of 
the rate of mass captured at apastron is needed to ob- 
tain such a transition. This is in principle feasible if an 
increment of the rate of mass lost by the Be star is ac- 
companied by an increase of the radius of the equatorial 
disc, so to exceed the maximum orbital separation (see 
the black dashed line in Fig. [1} . In this f ramework, the 
super-orbital variability hinted at in TeV ([Li et al.ll2012|) 
would be understood in terms of a transition towards the 
propeller state at all orbital phases, due to an increase 
of the mass loss rate of the Be star and by an increase of 
its size. Such a significant increase of the mass capture 
rate when the NS is close to apastron would be also able 
to explain the observed 7-ray luminosity in terms of the 
propeller emission [see Eq. (f2"9")) ]. 

7. DISCUSSION 

According to th e conventional picture (see, e.g. 
ILipunov et all 119921 ). a NS in a binary system evolves 
through ejector, propeller and accretor states, as it spins 
down. While on long timescales the transition between 
these states is mainly set by the evolution of the spin pe- 
riod of the NS, on shorter timescales a key role is played 
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by the rate at which the compact object captures the 
mass lost by the companion star. In particular, the large 
values of the ratio between the maximum and minimum 
rate of mass captured by the NS achieved if the orbit is 
eccentric, may induce state transitions along an orbital 
cycle, a flip-flop state. The range of mass capture rates 
spanned through an orbital cycle is even larger if the non 
degenerate star is of Be class, losing mass also through 
a dense equatorial disc which is transversed by the NS 
when it is close to pcriastron. The existence of systems 
alternating states on the timescale set by the orbital pe- 
riod is therefore a natural consequence of the way a NS 
evolves. 

While the observation of two magnetar-like bursts from 
a few-arcmin region compatible with the position of 
LS I +61°303 provided a strong indication in favor of 
a NS nature of the compact object in the system, it is 
unlikely that an accreting NS is hosted by LS I +61°303. 
Accreting NS in Be-HMXB show in fact X-ray pulsations 
at a period larger than few seconds. Moreover, their X- 
ray energy continuum spectrum is described by a power 
law with an exponential cutoff at 10-30 kcV, on which cy- 
clotron and / or iron emis sion features arc generally super- 
imposed (see IRei el 120 111 and references therein). No X- 
ray pulsation has been detected so far from LS I +61°303 
(|Rea et al.ll2010bD . while its X-ray spectrum is feature- 
less and does not show a cut-off in the X-ray energy band 
(jChernvakova et al.|[200l IZhang et al.ll2010t) . 

On the other hand, the luminosity observed from the 
system is of the same order of the spin-down power lib- 
erated by the NS when its spin period is close to the 
critical value marking the transition from the ejector 
to the flip-flop phase. An important clue to estimate 
the likelihood of finding a NS in LS I +61°303 in such 
a phase follows from the assessment of the time spent 
by the system in such a state. We showed how such 
a timescale crucially depends on the details of the as- 
sumed propeller torque. By assuming the rotating NS to 
release to the incoming matter the energy needed to un- 
bind it (see Eq. [T3|) . or the energy effectively stored in the 
NS rotation (Eq. IT4"]) . largely different estimates of the 
evolutionary timescales are obtained. This is essentially 
because the rotational energy of the NS when spinning 
close to the ejector/propeller transition largely exceeds 
the gravitational energy possessed by the in-falling mat- 
ter, evaluated at the large magnetospheric radius implied 
by a strong magnetic field. Among the propeller torques 
considered here, the weaker yields a total duration of the 
flip-flop phase which is much larger (sa 280 kyr) than the 
one implied by the stronger torque (« 25 kyr). Even con- 
sidering the stronger propeller torque, which is favored 
if the system is powered by the propeller effect when the 
NS is close to pcriastron, the expected total duration of 
the flip-flop phase [see Eq. (|27|) ] is larger by a factor of 
w 4 with respect to the timescale spent by the object in 
the pure ejector state: 



0.24 / • max\0.27 



(rh? ax f 



(30) 



where we have made explicit only the dependence of 
Eq. (HU and (J2ZJ) on the NS magnetic field, and on the 
maximum mass accretion rate. It is then reasonable to 
find the system in the flip-flop state. Moreover, even if it 



TABLE 4 

Ante-diluvian systems and LS I +61°303 



Name 



Ps(b) 



orb 



(d) 



Bi(G) < a > M 2 (M ) 



J 1740-3052 
J1638-4725 
J0045-7319 
B1259-63 

LS I +61°303 



0.57 
0.76 
0.93 
0.048 



231.0 
1941 
51.1 
1237 

26.5 



0.57 
0.95 
0.81 
0.87 

0.63 



3.9xl0 12 
1.9xl0 12 
2.1xl0 12 
3.3x10" 



11.0-15.8 

5.8- 8.1 

3.9- 5.3 
3.1-4.1 
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IManchester et all I I2U05 

( a ) The magnetic field of pulsars is determined as 3.2 X 10 19 (PP) 1 / 2 
G. 

is taken into account that the system must be relatively 
young to emit a spin-down luminosity £ few x 10 35 erg 
s _1 , the time spent by the system in the flip-flop phase 
is of few kyr, comparable to the total duration of the 
previous ejector phase, and yielding an age of the system 
f=a 5-10 kyr, of the order of typical spin-down ages of 
magnctars. On the other hand, a spin-down luminosity 
of the order of « 10 37 erg s _1 , would imply a smaller 
age for the system, with the NS emitting as an ejector 
all along the orbit. 

Either it lying in the ejector or in the flip-flop state, the 
presence of a young NS in LS I +61°303 would make the 
system closely related to so-called ante-diluvian systems 
(|van den Heuvelll200i[ ). rotation powered pulsars orbit- 
ing a high mass companion in a eccentric orbit, consid- 
ered as the progenitors of HMXB. The properties of the 
four sources of this class discovered so far are listed in 
Table H together with those of LS I +61°303. We note 
that B1259-63 is also one of the brightest 7-ray binaries 
known. All these NS have a superficial magnetic field 
in the range 3 x 10 11 -4 x 10 12 G, as derived from the 
observed spin down rate. Indeed, it was early proposed 
that some of these systems could be found in a propeller 
state when the NS was close to periastron. However, 
the rate at which the companion star would have to lose 
mass in order to overcome the pulsar pressure should be 
much larger than expected dCampana et aT1ll995l : iGhoshl 
[19951 iTavani fc Aronsl 119971 and references therein) or 
observed (jKaspi et al.lll996b[ ). Despite a proper evalua- 
tion of the ratio defined by Eq. ()30[) for the known ante- 
diluvian systems is far from the scope of this paper, it 
can be noted how the likelihood of observing a system in 
the pure ejector state increases when the maximum mass 
capture rate decreases, as it is the case of at least three 
out of the four known ante-diluvian systems, lying in a 
much larger orbit with respect to that of LS I +61°303. 

If confirmed by the discovery of pulsations, 
LS I +61°303 would be the first magnetar to be 
discovered in a binary system. The large luminosity 
variation shown by super fast X-ray transients on 
short timescales led IBozzo et alJ ()2008l ) to argue how 
a magnetar-like magnetic field could represent an 
efficient gating mechanism to make these systems to 
rapidly switch between propeller and accreting states. 
That a number of accreting HMXB should have host 
in the past a NS with a field in the magnctar-rangc 
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has also been claimed on the basis of their very large 
spin periods (e.g. the cases of 2S 0114+65 with a 
period of 2.7 hr, iLi k van den Heuvell 119991 and 4U 



2206+54 with a period of 55 6 s. iFinger et all l2010t 



llkhsanov k Beskrovnaval 120101 iReig et all l2012f ). The 
spin period at which a system eventually starts accreting 
mass increases in fact with the strength of the magnetic 
field (see Eq. Il6l an d the similar expression derived by 
IShakura et al.1 120121 for the equilibrium period of NS in 
the settling regime) , essentially because the value of the 
magnetic field sets the strength of the propeller torques 
that are responsible for the NS spin down. In this 
context the case of the Be/X-ray binary in the Small 
Magellanic Cloud, SXP 1062, with a period of 1062 s 
and with an estimated age of 16 ky r (iHaberl et al.ll2012ft 
is noteworthy. iPopov k Turollal ()2012h showed how 
such a short age strongly points to the presence of a NS 
with a large initial magnetic field in the system, ~ 10 14 
G, which could have spun the NS down very efficiently 
before the start of the accretion phase. A comparison 
of the age proposed for SXP 1062 with the timescales 
of the different propeller mechanism listed in Table [31 
indicates how the faster and stronger expression for 
the propeller torque is probably closer to the torque 
effectively experienced by a NS in that system during 
the propeller stage. How a strong propeller torque, of 
the order of that given by Eq. (fTlj) . should be possibly 
preferred in describing the evolution of systems like 



LS I +61°303, is also indicated by the ratio between the 
ejector and flip-flop timescales defined by Eq. ((301) . If 
the weaker propeller torques were in place, a NS would 
spend a much larger time in the flip-flop state than in 
a pure ejector state, and this is not indicated by the 
number of systems observed in the latter state (4) with 
respect to the single possible case of LS I +61°303. 

We finally note that, contrary to the case of SXP 1062, 
no association with a su pernova remnant could be been 
made for LS I +61° 303 (jFrail et aLlfl987h . This is not 
entirely surprising for a source with an estimated age 
between 10 and 20 kyr, since such an association can 
be found only for a fraction 0.64 and 0.55 of the radio 
pulsars with such ages estimated from their electromag- 
netic spin down, respe ctively (ATNF pulsar database, 
Manchester et al. 2005). Furthermore, no other Be X-ray 
binary has been observed embedded in a SNR, despite 
the relatively young age and large number of observed 
systems. This is most probably due to the large wind of 
the two progenitor massive stars, which has swept away 
most of the material around the binary, resulting in an 
under-luminous (hence undetectable) SNR after the ex- 
plosion. 

This work was supported by the grants AYA2009- 
07391 and SGR2009-811, as well as the Formosa pro- 
gram TW2010005 and iLINK program 2011-0303. We 
are grateful to Rosalba Perna and Jose Pons for useful 
discussions. 
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